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AERONAUTIC SYMBOLS 
I. FUNDAMENTAL AND DERIVED UNITS 



Symbol 

Metric 

English 

Unit 

Abbrevia- 

tion 

Unit 

Abbrevia- 

tion 

Length _ _ 

l 

meter . 

m 

foot (or mile) 

ft (or mi) 

Time _ 

t 

second 

s 

second (or hour) 

spc for hr') 

Force... 

F 

weight of 1 kilogram 

kg 

weight of 1 pound 

Ovu in J 

lb 

Power. 

P 

horsepower (metric) 


horsepower 

ho 


V 

/kilometers per hour. __ 

kph 

miles per hour _ _ 

mph 

Speed 

V 

\ meters per second 

mps 

feet per second 

fps 


2. GENERAL SYMBOLS 


Weight =ra <7 

Standard acceleration of gravity = 9. 80665 m/s 2 
or 32.1740 ft/sec 2 
W 

Mass= — 

9 

Moment of inertia =mP. (Indicate axis of 
radius of gyration k by proper subscript.) 
Coefficient of viscosity 


v Kinematic viscosity 

p Density (mass per unit volume) 

Standard density of dry air, 0.12497 kg-m~ 4 -s 2 at 15° C 
and 760 mm; or 0.002378 lb-ft“ 4 sec 2 
Specific weight of “standard” air, 1.2255 kg/m 8 or 
0.07651 lb/cu ft 


3. AERODYNAMIC SYMBOLS 


Area 

Area of wing 
Gap 
Span 
Chord 

Aspect ratio, ^ 

True air speed 
Dynamic pressure, ^pP 2 


Lift, absolute coefficient C L = 


<LS 

Drag, absolute coefficient 0©=-2 

qb 

Profile drag, absolute coefficient <7^=^ 
Induced drag, absolute coefficient C D4 =^ 
Parasite drag, absolute coefficient C Dp = 
Cross-wind force, absolute coefficient C c = 


C 

qS 


i* 

it 

Q 

12 

R 


a 

e 

(*0 

Vi 

Oi a 

y 


Angle of setting of wings (relative to thrust line) 
Angle of stabilizer setting (relative to thrust 
line) 

Resultant moment 
Resultant angular velocity 


Reynolds number, p — where l is a linear dimen- 

sion (e.g., for an airfoil of 1 .0 ft chord, 100 mph, 
standard pressure at 15° C, the corresponding 
Reynolds number is 935,400; or for an airfoil 
of 1.0 m chord, 100 mps, the corresponding 
Reynolds number is 6,865,000) 

Angle of attack 

Angle of downwash 

Angle of attack, infinite aspect ratio 

Angle of attack, induced 

Angle of attack, absolute (measured from zero- 
lift position) 

Flight-path angle 
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APPLICATION OF THE ANALOGY BETWEEN WATER FLOW WITH A FREE SURFACE 
AND TWO-DIMENSIONAL COMPRESSIBLE GAS FLOW 


By W. James Orlin, Norman J. Lindner, 
and Jack G. Bitterly 


SUMMARY 

The theory of the hydraulic analogy — that is, the analogy 
between water flow with a free surface and two-dimensional 
compressible gas flow — and the limitations and conditions o f the 
analogy are discussed. A test was run using the hydraulic 
analogy as applied to the flow about circular cylinders o f various 
diameters at subsonic velocities extending into the supercritical 
range. The apparatus and techniques used in this application 
are described and criticized. Reasonably satisfactory agree- 
ment of pressure distributions and flow fields existed between 
water and air flow about corresponding bodies. This agreement 
indicated the possibility of extending experimental compressi- 
bility research by new methods. 

INTRODUCTION 

An analogy exists between water flow with a free surface 
and two-dimensional compressible gas flow (hydraulic 
analogy). The water must flow over a smooth horizontal 
surface bounded by vertical walls geometrically similar to 
the walls bounding the corresponding compressible gas flow. 

The mathematical basis of this hydraulic analogy was 
presented by Riaboucliinsky in reference 1, in which he also 
described his apparatus for investigating the flow in a Laval 
nozzle. In reference 2, he extended the theory to include 
drag considerations and outlined the probable usefulness of 
the hydraulic analogy. Binnie and Hooker, in reference 3, 
obtained surveys along the center line of a channel with a 
constriction. By employing the characteristics method to 
calculate accurately the flow in a Laval nozzle, Preiswerk, 
in reference 4, demonstrated conclusively that the methods 
of gas dynamics can be applied to water flow with a free 
surface. 

The National Advisory Committee for Aeronautics became 
interested in the hydraulic analogy because it seemed an 
easy and inexpensive way of studying two-dimensional com- 
pressible gas flow; in particular, phenomena occurring in air 
at speeds too high for visual observations could be observed 
at very low speeds (3 or 4 fps) in a water channel. Pre- 
liminary investigations were made in the Langley tank no. 1 , 
where difficulty with the vertical accelerations — assumed 
negligibly small in the analogy — was experienced. 

A water channel was designed and constructed in the 
Langley 8-foot high-speed tunnel in the spring of 1940. The 


channel was so constructed that flow fields involving both 
subsonic and supersonic velocities about aerodynamic bodies 
could be investigated. The value of the analogy in such 
flow fields has not been previously demonstrated. The 
development of the measuring apparatus and techniques is 
presented herein. The application of the analogy to flows 
through nozzles and about circular cylinders at subsonic 
velocities extending into the supercritical range is also 


presented. 

SYMBOLS 

Cp 

specific heat at constant pressure 

Cv 

specific heat at constant volume 

7 

adiabatic gas constant, ratio of c v to c v of gas 

T 

absolute temperature of gas 

p 

mass density of gas 

P 

pressure of gas 

h 

enthalpy or total heat content 

<Z 

dynamic pressure of gas pV 

<j 

surface tension of liquid 

P 

pressure coefficient 

p 

viscosity of liquid 

/ i — \ 

a 

speed of sound in gas 7 

V 

velocity of flow 

d 

water depth 

M 

Mach number, stream value unless otherwise 
indicated ( V/a for gas; V/fgd for water) 

9 

acceleration of gravity 

F c 

compressibility factor 

x, y 

rectangular coordinate axes 

U, V 

components of velocity in ^-direction and 
^/-direction, respectively 

<t> 

velocity potential in two-dimensional flow 

X 

wave length of surface waves in fluid 

U 

velocity of propagation of surface waves in 
fluid 

D 

diameter of circular cylinder 

R 

Reynolds number ( VDp/p ) 

e 

angular measurement clockwise about cylinder; 
at stagnation point, 6= 0° 


l 
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Subscripts: 

No subscript any value of variable 

0 value at stagnation (F=0) 

st value for completely undisturbed static con- 

ditions, no flow in channel 

1 local value of variable; value at surface of 

model, at channel walls, or in field of flow 
s value in undisturbed stream 

1, 2 any two values of variable 

max maximum value of variable 

cr critical value of variable, value at which Mach 

number of 1 is reached at some point in flow 
x partial derivative with respect to x 

(for example: 4>,=^; 

y partial derivative with respect to y 

THEORY 

The following is a condensation of the theory and mathe- 
matical development of the hydraulic analogy as given by 
Preiswerk (reference 4 ). 

Two assumptions are made in the mathematical develop- 
ment: 

(1) The flow is irrotational. 

(2) The vertical accelerations at the free surface are 
negligible compared with the acceleration of gravity. The 
pressure in the fluid at any point therefore depends only on 
the height of the free surface above that point. 

The analogy between the flow of water with free surface 
and the flow of a compressible gas may be obtained by set- 
ting up the energy equations for each. From the energy 
equation for water, the velocity is 


From equations ( 4 ) and ( 5 ), a further condition for the 
analogy may be derived 

d =P 

do Po 


Since for adiabatic isentropic flow in the gas 

i 

P = 

Po 


TV-i 

To) 


and since, in the analogy, from equations ( 3 ) and (6) 


p = rf = T 
Po do T 0 


then 


T / T 


i 


To \To) 

and, therefore, the analogy requires that 7= 2.0 
From the relation 

y 

Po \Po/ 

it may be seen, since 7 = 2.0, that 

Hi) 


Po 


(6) 


17 ') 


( 8 ) 


The velocity potential for water is given by the equation 
4,11 (' ~%) + 4>! " ; 41 (id 


***■•'= 0 ( 9 ) 


and 


V*= 2 g{d 0 -d) 


1 max A “ f/d(l 


(la) 

(lb) 


and the corresponding equation for a gas is 


0XX ^1—^2 ) + 0i/«/ (l 


0 /01/ 


= 0 


( 10 ) 


The corresponding equations for the velocity of a gas are 


Therefore, for identical expressions, 


and 


V 2 = 2 g(h 0 —h)= 2 gc p (T 0 — T) 
1 max \ 2 g/i 0 \[ 2 gc p T 0 


(21 


Therefore, if the ratio V/V m ax for g&s I s equated to V/ V 7h 
for water, 

do d K h 1 0 1 

do ho To 


or 


d = T 

do 1 0 

The equation of continuity for water is 

d(w/) d(wO 
dx dy 

The continuity equation for two-dimensional gas flow is 
d( 7 /p) d(op) =() 


dx 


dy 


gd __ a 2 
2 ga 0 2 gh () 


( 11 ) 


From equations ( 9 ) and (10), the velocity ^gd in (he 
liquid flow is seen to correspond to the velocity of sound in 
gas flow. In a subsequent section entitled “Discussion,” the 
value of ^fgd is shown to be the velocity of propagation of 
surface waves, the wave lengths of which are large in com- 
parison with the water depth. The ratio of 1 / -yfgd in the 
liquid flow corresponds to the Mach number Vja in the 
gas flow. 

If the velocity of the liquid flow is less than yW (M<M), 
the water is said to be “streaming.” If the velocity of 
liquid flow is greater than \gd the water is said 

to be “shooting.” 

In shooting water under certain conditions, the velocity 
of the flow may strongly decrease for short distances and 
the depth may increase. A11 unsteady motion of this type 
is called a hydraulic jump. Hydraulic jumps of small 
intensity are propagated with the velocity ifgd. 
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The analogy may be summarized as follows: 


Significant quantities and characteristics 
of two-dimensional compressible gas 
flow, 7=2 

Corresponding values in analogous liquid 
flow 

Temperature ratio, 77 To 

Density ratio, p/p 0 . 

Pressure ratio, pf po - 

Velocity of sound, a = -/ 7 P 

V p 

Mach number, V/a 

Subsonic flow 

Supersonic flow 

Shock wave 

Water-depth ratio, d/d 0 
Water-depth ratio, d/do 
Square of water-depth ratio, (d/do) 2 

Wave velocity, -yjqd 

Mach number, V/^gd 
Streaming water 
Shoot ing water 
Hydraulic jump 


The Mach numbers in the liquid flow may easily be com- 
puted. If equation (la) is substituted in the expression for 
stream Mach number 


M= 


Vs 


the expression becomes 


(12a) 


The points for M—l are located where the depth is two- 
thirds of the total head, as is easily verified from equa- 
tion (12a). The local Mach numbers are computed from 


M t = 


F2 (do— r/ ( ) 

L d, 


] 


1/2 


(12b) 


The pressure coefficient, a significant quantity in com- 
pressible gas flow, may be computed for liquid flow. If 


m - 


is substituted for in the expression for pressure 
V* 

coefficient 


HP 


the expression becomes 


P=F t 


‘(DH 

.(£)'- *J 


For a compressible gas 


(13) 




7 


(i+ 7 -Hm/) 7 "'-i 


and with 7 = 2.0, as required by the hydraulic analogy, the 
compressibility factor is exactly 

f c =i+\m? 

When the analogy is applied to the study of air flow, 
accurate quantitative results will not be obtained because, 
for strict agreement between water flow and gas flow, 7 
must equal 2.0, whereas for air 7 is 1.4. The relation 
between pressure ratio and local Mach number for the two 


values of 7 is shown in figure 1. A given pressure ratio 
corresponds to a higher Mach number in air than in the 
fictitious gas with 7 = 2.0. The corresponding depth ratios 
d/d 0 shown in figure 1 are equal to the square root of the 
pressure ratios for 7=2.0. With a given stream Mach 
number, the value of the critical negative-pressure coefficient 
(fig. 2) is greater in air than in the gas with 7 = 2.0. The 
differences are not very large, however, and the flow phe- 
nomena observed in the water flow should be qualitatively the 
same as those occurring in the two-dimensional compressible 
flow of air. 

APPARATUS 

The tests were conducted in the vertical, return-flow 
water channel which was designed and constructed in the 



Figure 1 .— Variation with local Mach number of pressure ratio p/po in a comp ressible gas, 
and the water-depth ratio d/do in the hydraulic analogy corresponding to Vp/Po for 7=2.0. 
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Langley 8-foot high-speed tunnel. Figure 3 shows the flow 
circuit and plan view of the free water surface and also the 
orifice locations of the water channel. A motor-driven 
propeller of 1-foot diameter forces the water into the ex- 
panding section and through a screen into a large quieting 
section. The fluid is accelerated in the convergent section 
and flows through the test section downstream to the propel- 
ler so that continuous flow is maintained. The test section 
could be replaced by the nozzles shown in figures 4 (a) to 
4 (c). The shape of the channel test section is given for 
comparison in figure 4 (d). The entire floor of the entrance 


5'0" 



Model position 
Or/fices on £ / 


:Olf fuser 




k, 


life 


— * 

s 

J 1 






Wo// pressure orifices 


'Total-head orifice 


Ipl) 



» Screen 
40x40 
mesh 


Section A~A 



Section B-B; C-C 


Section D~D 


Section E~E\ E~E 


Screen 

80x80 

mesh 

0 / 


Scale, ft 


Section G~G 
( c ) 

(a) Plan view. 

(b) Elevation. 

(c) Cross sections. 



Figure .‘1. — Schematic view of water channel with sections. 


cone and test section is horizontal except for the diffuser 
which has a 1° slope. For a static depth of 1.5 inches, the 
maximum volume flow is approximately 540 cubic inches per 
second and the power required is 1/10 horsepower. 





(a) Two-inch Laval nozzle (reference 4). 

(h) Four-inch Laval nozzle (reference 4). 

(c) Four-inch Laval nozzle with modified approach. 

(d) Twenty-four-inch nozzle, final wall shape. 


Figure 4.— Nozzle shapes for depth-ratio flow tests. 


AN ANALOGY BETWEEN WATER FLOW AND GAS FLOW 


5 


Water depths at points in the test section were obtained 
with the survey equipment shown in figures 5 and 6. A 
block on which was mounted a vertical shaft tipped with a 
fine wire probe (0.030-inch diameter) could be moved both 
parallel and perpendicular to the channel axis. The probe 
was capable of 2%-inch vertical travel continuously measur- 
able by a standard micrometer to the nearest one-half 
thousandth of an inch. A small neon lamp was mounted 
in series with the insulated probe and the water to provide a 
positive sign of contact between the probe and water level. 

A similar micrometer unit, movable along a single fixed 
rail, was located upstream along the floor center line to 
measure total head. Water depths could also be measured 


by static orifices on the channel walls and along the floor. 
Each orifice was connected by a valve to a common sump 
and burette. By opening only one valve at a time, individual 
pressures could be obtained in the burette by means of a 
vernier height gage fitted with a hairline mirror sight. The 
locations of these orifices are shown in figure 3. 

Photographs and motion pictures of the free surface could 
be obtained by cameras mounted above the test section. A 
Strobolux and Strobotac were used as the source of illumina- 
tion. A circular glass plate, fitted into the center of the floor 
of the test section, permitted illumination of the water 
surfaces from below. 

METHODS 



A Longitudinal rail D Cross rail 

B Micrometer head E Longitudinal rail 

C Survey probe F Model 

G Glass plate 

Figure 5.— Test section and survey equipment seen from downstream end of channel. 



A Strobolux 
B Strobotac 
C Camera 
D Micrometer head 


E Survey probe 

F Support for ground-glass screen 
G Mirror 

H Position of light source 


The two-dimensional models were mounted at the center 
of the glass plate, one end flush with the plate and the other 
end extending above the water surface. Orifices were drilled 
in the model a fixed distance above the glass plate and normal 
to the model surface. 

From the theory of the hydraulic analogy, it is readily 
seen that accurate measurements must be made of the depths 
in alt measurable points in the channel. For each operating 
speed two readings were usually obtained for each datum 
point, one by the burette and the other by the survey 
carriage. The burette reading, by means of the hairline sight 
and vernier, could be taken as quickly or as often as neces- 
sary without affecting stream conditions; whereas, because 
of the formation of capillary and standing waves, the surface 
probe reading was good only at the instant the probe touched 
the water. This effect was not serious at local supersonic 
speeds inasmuch as the standing waves were then unable to 
move upstream from the point of water contact. When the 
disturbance was removed by breaking contact with the water 
surface, original stream conditions were restored in about 
10 seconds, at which time another reading could be made. 

The micrometer unit movable along a single fixed rail 
measured the depth at the upstream positions and the depths 
at the first pressure orifice on the channel center line. From 
continuity considerations, the total head 


d*=(l s 


(L — r/„ 



(14) 


where b is the breadth of the channel and subscript 1 refers 
to the station at the first orifice on the center line. The 
stream depth d s was measured sufficiently far ahead in the 
test section to be unaffected by the presence of the model. 
The total head, however, could be measured directly by the 
upstream micrometer in the vicinity of the boundary wall 
at the upstream end (fig. 3). The upstream micrometer was 
calibrated against the test-section survey probe before each 
run. All depths were measured from the level floor of the 
channel. The net accuracy of all readings was within 
±0.002 inch. 

Orifice pressures were used to obtain the water depth on 
the surface of the models and at the test-section side walls. 


Figure (>.— Approach and tost section showing equipment for photography. 
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This typo of measurement is necessary inasmuch as the 
capillary rise of water on all bodies that break the water 
surface invalidates the use of the survey probe. The probe 
is of great value, however, in obtaining surveys at points 
where the influence of capillarity is negligible, which in most 
cases is about ){ inch from an exposed surface. Extreme care 
had to be taken in all pressure measurements to hold capillary 
rise constant in the burette, inasmuch as very slight concen- 
trations of dirt on the glass would change the meniscus shape. 
The burette was, therefore, cleaned at regular intervals with 
chromic acid, and a 10-percent solution of aerosol was placed 
on the meniscus to reduce surface tension and capillarity. 
This procedure enabled consistent readings within the desired 
accuracy. 

The shadowgraph still photographs were obtained by a 
5- by 7-inch view camera. The test setup and schematic 
diagram are shown in figures 6 and 7 (a). In the initial 
trials it was found that the most convenient light source, 
with existing collimator lenses, was a Strobolux-Strobotac 
unit; and by adjusting th(‘ flash frequency, various light in- 
tensities could be obtained. The refraction patterns formed 
by a vortex, hydraulic jump, and capillary wave are shown 
in figures 7 (b) to 7 (d); and from these patterns or combina- 
tions, a qualitative analysis of shadowgraphs may be made. 
In order to reduce excessive pattern refraction, the ground- 
glass screen was placed as close to the water surface as pos- 
sible. Motion pictures of wake frequencies and alternate 
trailing vortices were obtained by the same methods with no 
great difficulty. 


L- Comer a lens 

h> 

I 

I 


. , - Grounc/-g/ass 

screen ... 


Cylinder mo de A. v | 
Channel floor 


Gloss plate-- I 


Parallel 
J light 


/Source 



Collimator 







7.— Analysis of method for photographing free surface. 


RESULTS 

In order to ascertain the nature and the applicability to 
air flow of data obtained in the water channel, tests were 
made with flow through nozzles and about circular cylinders. 
These tests were run at various depths and Mach numbers 
and with models of various sizes in the channel. The results 
are first presented without discussion. 

Plots of the water depths and water-depth ratios along the 
center line of convergent -divergent nozzles as measured by 
the test-section micrometer unit are shown in figures 8 to 10. 
The shapes of these nozzles are shown in figures 4 (a), 4 (b), 
and 4 (c), respectively. The curves labeled “Reynolds flow” 
are computed from one-dimensional theory; that is, the ve- 
locities at all points in a cross section are assumed the same 
in magnitude and parallel to the center line of the channel. 
The depth at any point is computed from the liquation of 
continuity 

T T X Area = Constant (15) 

where 

Area = Breadth of the cliannelXWater depth 

The water-depth ratios are seen to increase with the increase 
in depth in figure 8. The nozzle in figure 10 was modified 
to have a very gradual approach section and the water-depth 
ratios for this nozzle are the same at any depth. 



(a) Variation of water depth along center line showing shape of free surface. 

(b) Variation of water-depth ratio d/do along center line for different values of the total 

head do. 

Figure 8.— Effect of change of total head on flow through a convergent-divergent nozzle. 
Throat width, 2 inches. Run at choking Mach number. Refer to figure 4 (a) for shape. 


Figure 
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(a) Variation of water depth along center line showing shape of free surface. 

(b) Variation of water-depth ratio d/do along center line for different values of the total 

head do. 


Figure 9. — Effect of change of total head on flow through a convergent-divergent nozzle. 
Throat width, 4 inches. Run at choking Mach number. Refer to figure 4 (b) for shape. 


Tli o variation with total head of the water depths and 
water-depth ratios about a circular cylinder as measured by 
means of the pressure orifices and burettes is shown in figure 1 1 . 
In this case the depth ratios decrease with increase in water 
depth, A comparison between micrometer and burette 
measurements is shown for a 2-inch nozzle in figure 12, where 
it is seen that the micrometer readings are greater than the 
burette readings and the difference increases with depth. 

Figures 13 and 14 show the variation with water depth of 
the shadowgraphs of the flow about a 2%-inch-diameter 
circular cylinder at Mach numbers of 0.40 and 0.60, respec- 
tively. The corresponding variation of pressure at 0=90° is 
also shown. Add i tional shadowgraphs of the flow about circu- 
lar cylinders of various sizes and with various water depths are 
shown in figures 15 to 19. The pressure distributions about 
the 2%-inch circular cylinder in the water channel with a 
static depth of % inch are compared for various Mach numbers 
in figure 20 with the pressure distributions about a %-inch 
circular cylinder in air. The data for the circular cylinder 
in air were available from tests made in the Langley rectangu- 
lar high-speed tunnel. 

With various sizes of circular cylinders in the test section 
the variation of the local Mach numbers along the channel 
wall for different stream Mach numbers is shown in figures 21 
to 24, and the corresponding stream depths d s and total 
heads d 0 are given in figure 25. 

DISCUSSION 



Distance along nozzle center Hne } in. 



(a) Variation of water depth along center line showing shape of free surface. 

(b) Variation of water-depth ratio d/do along center line for different values of the total 

head. do. 


Figure 10.— Effect of change of total head on flow through a convergent-divergent nozzle. 
Throat width, 4 inches (modified approach section). Run at choking Mach number. 
Refer to figure 4 (c) for shape. 


channel, and the effect of these factors must be ascertained 
before the channel can be used for the investigation of the 
analogous gas flows. 

From an examination of figures 15 to 19 it is seen that 
concentric waves appear in front of the model. The fact 
that these waves do not appear in the schlieren photographs 
(fig. 26) of air flow taken at the Langley rectangular high- 
speed tunnel suggests that these waves are not part of the 
analogy. 

A disturbance at some point in a liquid generally will give 
rise to two types of waves — short surface-tension waves 
called capillary waves and considerably longer gravity waves 
(reference 5 and pp. 353-402 of reference 6). The expression 
for the velocity of propagation of water waves is 


U= 



+ 


2tt(A 
PA / 


tanh 


27rdJ /2 


(16) 


Figure 27 is a plot of V against A for water depths of 0.5 
inch and 2.0 inches. This figure shows that waves cannot 
exist with a velocity of propagation less than approximately 
0.75 foot per second. Disturbances of wave length less than 
that corresponding to the minimum wave velocity are termed 
capillary waves since they depend primarily on the surface 
tension of the fluid. Propagation velocities greater than 
this minimum velocity correspond to shorter capillary waves 
and longer gravity waves. 

If (7=0 and g?<C A, equation (16) becomes 


Various factors influence the results obtained in the water 


U=Jgd 


(17) 


823737—49 2 
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O 10 20 30 40 50 60 70 80 90 


6, deg 

(a) Water depths as determined by burette. 

(b) Variation of water-depth ratio rf/rfo as determined by burette for different values of 

total head. 

Figure li.- -Effect of change of total head on burette measurements at the surface of a 
2!4-inch-diameter circular cylinder. A/=0.70. 


and t his velocity — the velocity of propagation of long grav- 
ity waves— is the basic surface wave velocity. The other 
waves defined by equation (16) are water surface waves 
and are not considered in the hydraulic analogy. Figure 
27 shows that changes in the water depth have little 
effect on the velocity of propagation of waves for wave lengths 
less than 0.1 foot. 

The wave lengths of some of the standing waves appearing 
in front of the models in figures 15 to 19 were measured from 
the photographs and were plotted in t lit' inset of figure 27 
against the local velocity (equal, for standing waves, to the 
velocity of propagation) at the point where the wave occurred. 
This plot showed that these waves were capillary waves. 
The check marks on the inset show the number of readings 
that fall on the plotted points. These capillary waves, 
although they have no part in the hydraulic analogy, have 
two adverse effects on the application of the analogy, namely, 



(a) Variation with total head do of water-depth ratio measured by burette and 

micrometer. 

(b) Comparison of micrometer and burette readings of depth. 


Figure 12.— Comparison of micrometer and burette measurements at the minimum section 

of th.* 2-inch nozzle. 


complication of the flow photographs and decrease in the 
accuracy of the depth measurements taken with the test- 
section probe. 

The results obtained in the water channel depend to a 
large extent on the depth at which the tests are run. The 
effect of depth on the formation of standing capillary waves 
is shown in figure 28. If high stream Mach numbers without 
capillary waves in the stream are to be obtained, a shallow 
depth is prescribed. 

The effect of depth on the flow patterns and pressure coeffi- 
cients about circular cylinders is very pronounced (figs. 13 
and 14). The pressure coefficient rises to a maximum value 
at r/ s ^0.9 inch. The decrease in pressure coefficient at 
lower depths is believed to be due to the effect of bottom 
boundary layer; the decrease at higher depths is due to the 
effect of the vertical accelerations. Comparison of figures 13 
and 14 with figure 26 indicates that the water flow most 
nearly corresponds to air flow at depths that are between 
0.75 inch and 1.0 inch. 
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(a) d,t= 0.50 inch. 


(d) d«t=1.50 inches. 



(a) d, t = 0.50 inch. 


(d) d» t = 1.50 inches. 




(b) d,i= 0.75 inch. 


(e) d*t=2.00 inches. 


(b) d,f=0.75 inch. 


(e) rf„<=2.00 inches. 



Figure 13.— Variation with static depth of flow pattern and pressure coefficient for the 
2 |4-inch -diameter circular cylinder. M=0.40. 


Figure 14.— Variation with static depth of flow pattern and pressure coefficient for the 
2}4-inch*diameter circular cylinder. M=0.60. 
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(.0 A/— 0.23. 


(d) A/=o.74. 



(a) A/=0.41. 



(d) A/ =0.00. 





(h) A/=0.50. 


(e) A/=G.7l. 


(e) Ai =0.02. 




(C) M=0.55. 


(f) A/=0.77. 


(f) A/=0.98. 




Figure 15.— Shadowgraphs of flow about the 0.30-inch-diameter circular cylinder. 
Static water depth, 1.00 inch. 


Figure 16.— Shadowgraphs of flow about the 1 ^-inch-diameter circular cylinder 
Static water depth, 0.76 ineh. 
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(a) M= 0.34. 


(d) A/=0.f>0. 


(a) A/=0.41. 


(d) Af=0.Gl. 



(b) Af= 0.51. 


(e) A/=0.G5. 



(b) Af=0.46. 


(0) A/=0.72. 




(c) Af=0.50. 


CO A/=0.78. 


(c) Af=0.55. 


(0 Af=0.71. 


Figure 17.— Shadowgraphs of flow aboutj'thc F/j-inch -diameter circular cylinder. 
Static water depth, 1.00 inch. 


Figure 18.— Shadowgraphs of flow about the 2U -inch-diameter circular cylinder. 
Static water depth, 0.75 inch. 
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(c) Af=0.4<J. (f) M=0.71. 

Figure 19.— Shadowgraphs of flow about the 2)4 -inch-diameter circular cylinder. 
Static water depth, 1.00 inch. 


The water is accelerated from the total head where the 
fluid velocity is zero to a value of zero acceleration and maxi- 
mum velocity in the test section. The motion of the water 
through this cycle is actuated by the force of gravity which 
is thus the important factor in the hydraulic analogy. If an 
obstruction is placed in the path of a fluid at constant veloc- 
ity, local accelerations of flow must take place in the field 
about the obstacle. If these accelerations are large, the 
vertical components, which have no part in the analogy, are 
no longer negligible in comparison with the acceleration of 
gravity; and a distortion of the flow results. This distortion 
is such that in a region of increasing deceleration, such as 
the region in which the (low is approaching stagnation, the 
depth is usually less than that required in the analogy; and 
the depth indicated by the pressure in the static orifices is 
greater than the true depth. In a region of increasing 
acceleration the opposite effects occur. The vertical accel- 
eration depends upon the water depth and upon the slope 
and curvature of the free surface. 

If the water depth is increased for two cases of the same 


stream Mach number, the slope and curvature of the free 
surface are likewise greater. This variation is shown in the 
following table in which the value of d s is computed from 
equation (12a) for chosen values of M and r/ 0 : 


M 

do 

d t 

d./do 

* ’ 

l 1 
■e 

| 

\ ’ 

(in. /sec) 

0. 920 

1 . 000 

0. 700 

0. 700 

0. 300 

10. 78 

. 926 

1. 100 

. 770 

. 700 

. 330 

11.29 


The slopes and curvatures of the free surface are approxi- 
mately proportional to the term do — d s , and the higher 
value of that term therefore corresponds in general to greater 
vertical accelerations. An illustration of the foregoing 
statements is given by the results of the nozzle tests in which, 
as the water depths decrease, the depth ratios approach 
those calculated by Reynolds flow (fig. 8); all tests were 
made at choking Mach number. Binnie and Hooker 
observed the same effect (reference 3). If the slope of the 
free surface is small and if the nozzle is so shaped that the 
acceleration in the entrance can be kept sufficiently small, 
no variation of water-depth ratios occurs with changing 
depth (fig. 10). The same effect is shown in figures 1 1 and 
12 as measured by static orifices and burettes. 

The effect of vertical accelerations on depth ratios deter- 
mined by means of static orifices is opposite to that obtained 
by the probe. It may be assumed therefore that this ad- 
verse effect would be negligible when the two methods of 
reading pressures coincide. The readings taken at the mini- 
mum section of the 2-inch nozzle (fig. 12) coincided when the 
depth was reduced to 1.3 inches. If the free surface slope 
had been more abrupt (smaller model or model having a 
sharp pressure rise), a still shallower depth would have been 
required to attain this condition; therefore, large models 
and small depths should be employed. The minimum depth 
is limited, however, by the effects of the boundary layer on 
the floor of the channel. 

The development of the boundary layer produces a velocity 
gradient along the channel similar to the gradient in a wind 
tunnel. This effect can be compensated in the same way by 
diverging the walls or in the case of the water channel by 
sloping the floor downward in the direction of flow. An effect 
for which compensation does not appear feasible is the 
thickening and thinning of the boundary layer in regions of 
decelerating and accelerating flow about a model, behavior 
which causes a distortion of the velocity distribution such as 
would correspond qualitatively to a thinning of the model. 
These boundary-layer effects are particularly serious in the 
water channel because of the appreciable boundary-layer dis- 
placement thickness (found to be of the order of 0.1 in.) in 
comparison with the total depth of the fluid. The boundary- 
layer effects could obviously be minimized by increasing the 
water depths, both because of the increased ratio of total 
depth to boundary-layer displacement thickness and because 
of the increased Reynolds numbers corresponding to the 
higher velocities with given Mach numbers. This require- 
ment is, unfortunately, inconsistent with the minimizing of 
the effects of the capillary waves and of the vertical 
accelerations. 



Pressure coefficient , P Pressure coefficient , P 
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(a) M =0.41. (b) Af =0.50. 

(c) A/=0.61. (d) M= 0.71. 

Figure 20.— Comparison of pressure distributions about a 2J4-inch circular cylinder in the water channel with the pressure distributions about a H’-ineh circular cylinder in the 

Langley rectangular high-speed tunnel (RUST). 
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Figure 21. -Variation of local Mach number along channel wall with a 0.3-ineh-diamcter 

circular cylinder. 


The analogy is further limited by the fact that it applies 
for a hypothetical gas having a value of 7 = 2 . 0 ; the effect of 
this limitation has not been completely determined. The 
effects shown in figures 1 and 2 are not very large, however, 
and can be taken into account in the interpretation of data 
taken in the water channel. Other investigations suggest that 
tin' influence of the value of y in subsonic compressible flow 
is not great. Kaplan (reference 7) thus found that to the 
third approximation the effect of the value of y in the com- 
pressible flow up to the critical speed was negligible and 
Von Ivarman (reference 8 ) obtained a widely used expression 
for change of pressure coefficient with Mach number by using 
the assumption 7 = — 1 . For supersonic flow, the character- 
istics curves are considerably influenced by the value of 7 . 
The ratio of maximum velocity to the velocity of sound is 
greater with 7 = 1.4 than with 7 = 2.0; and for a given change 
in flow angle, the velocity change, as determined by the 
characteristics method, is also greater. 

The Reynolds numbers in the water channel are calculated 
for complete submersion and are very low; the maximum 
value of R for the 2 ^-inch-diameter cylinder with a static 
water depth of 0.75 inch is 15,000. The comparison of the 
water-channel and wind-tunnel data (figs. 15 to 19 and 26) 
showed that, although the cylinder tests were made at 



Figure 22. — Variation of local Mach number along channel wall with a 1 ^-inch-diamctcr 

circular cylinder. 

considerably higher Reynolds numbers in air than in water, the 
same type of flow occurred in both cases; the main consid- 
eration is agreement of the flow patterns rather than identical 
Reynolds numbers. 

In a wind tunnel, a choking or maximum stream Mach 
number occurs when local Mach numbers of 1.00 extend 
across the section between the model and tunnel walls. A 
similar effect was noted in the water channel. Figure 25 
shows that the' stream depth decreases and the total head 
increases with Mach number up to a maximum value that 
depends on the size of the model tested. At this choking 
condition, both the stream deptli and the total head increase 
with an additional power input. Figure 29 shows the maxi- 
mum Mach number for various ratios of cylinder diameter to 
channel width. Variation of the data from the theoretical 
curves might be expected inasmuch as the theoretical values 
were computed by one-dimensional theory (Reynolds flow); 
whereas, the actual flow is two-dimensional. Maximum 
Mach numbers that occur in wind tunnels, however, agree 
much more closely with the one-dimensional theory than do 
those shown in figure 29 (see reference 9); and it therefore 
seems likely that the divergence is due to other factors, such 
as boundary-layer effects and vertical accelerations, rather 
than to the two-dimensional nature of the flow. 
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Figure 23.— Variation of loctil Mach number along channel wall with a 2}4 -inch-diameter 

circular cylinder. 



Figure 24.— Variation of local Mach number along channel wall with a 6-inch-diameter 

circular cylinder. 


The disturbance along the channel wall with various size 
cylinders in the test section is shown in figures 21 to 24. 
The wall disturbance is an indication of the severity at 
high Mach numbers of channel-wall or wind-tunnel-wall 
interference. As the stream Mach number approaches 1.00, 
the wall disturbance increases until at the choking condition 
the interference becomes large even for the smallest cylinder. 
For the 1 K-inch-diameter cylinder (fig. 22), the point of 
sonic velocity at the wall is located approximately 3 diameters 
behind the cylinder. This result is in agreement with the 
Mow photographs (figs. 15 and 17), which show a gull-wing- 
shaped shock starting from the region behind the cylinder 
and extending downstream toward the walls. Figure 30 
shows the lines of constant Mach number in the Mow field 
about a 6-inch-diameter cylinder at choking Mach number 
and shows clearly the supersonic Mow occurring behind 
the cylinder. 

A comparison of the shadowgraphs (figs. 15 to 19) with 
schlieren photographs of Mow about circular cylinders in air 
(fig. 26) shows that the two flows are very similar. The type 
of Mow is the same in the two cases and is the type that is 
characteristic of Reynolds numbers somewhat below the 
critical value. At Reynolds numbers between 50 and 
350,000, an unstable condition is set up behind a circular 


cylinder in which vortices are shed alternately from each 
side. Laminar separation occurs at a point approximately 
80° from the forward stagnation point and a vortex sheet is 
formed which extends downstream and finally rolls up into 
a large vortex. The vortices so shed arrange themselves 
into a Karman street. (See reference 6, pp. 217-218.) 



Figure 25. — Variation with stream Mach number of stream depth and total head with 
circular cylinders in the test section. rf,*=l.()0 inch. 


Velocity of propagation \ U ) fps 


16 


REPORT NO. 875 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



(a) Af=0.401; /2 = 21<;,00(); IJ= 1.0 inch. (d) A/=0.730; 22=337,000; Z)=1.0 inch. 



(b) Af= 0.600; 22 = 290,000; D= 1.0 inch. (e) A2=0.880; = 73,0(X); I) =0.2 inch. 



(c) At =0.629; 22 = 308,000; 2)=1.0 inch. (f) A/=0.880; 22 = 73,000; D= 0.2 inch. 

(Same as (e) showing extended field.) 


Figure 26.— Schlieren photographs of circular cylinders in the Langley rectangular high- 
speed tunnel. 



Figures 16 (a), 18 (a), 19 (b), and 26 (a) show this type ol 
flow, which is the same in the water channel as in the air 


flow. Because of the occurrence of separation, the actual 
pressure distribution is different from that calculated by 
potential-flow theory. The actual pressures over the for- 
ward part of the cylinder are higher and those over the rear 
are lower than those calculated by theory. The negative 
pressure peak is greatly reduced. Figure 20 shows the pressure 
distributions about circular cylinders in air and in the 
water channel. The quantitative results obtained in the 
water channel are very close to those obtained in air 
(figs. 20 (a) to 20 (c)). This close agreement is believed, how- 
ever, to be largely fortuitous, resulting from an interaction 
between the bottom-boundary-layer effects and the rather 
large channel-wall interference in the water channel. At the 
choking Mach number M= 0.71 in the water channel (see 
fig. 20 (d)), the interference effects exceed the relieving 
effects of the thinning of the boundary layer, and the nega- 
tive pressure coefficients near 0=90° exceed those obtained 
in the wind tunnel. 

Because of the unsteady nature of the flow, instantaneous 
velocities greater than those indicated by the pressure dis- 
tributions occur near the surface of a cylinder. The critical 
speed may therefore be expected to be lower than indicated 
by the pressure distributions, an effect that has been observed 
in air. Figure 2 shows that with a given peak negative 
pressure coefficient the critical Mach number should be 
somewhat lower in the water channel (y = 2. 0) than in air 
(7=1.4). 

Another effect of the unsteady flow is the asymmetrical 
development of shock. When the flow closes in on one side 
of the cylinder, the induced velocities increase and the local 
Mach number may exceed 1.00. A shock wave may thus 
form on one side while none exists on the other or a more 
intense shock may exist on one side than on the other. Such 
asymmetrical shock patterns are shown in figures 16 (b) to 
16 (d), 17 (b), 18 (b) to 18 (d), and 26 (b) to 26 (d). Exactly 
the same phenomenon occurs in the water channel as in the 
air flow. In some cases, the shedding of these waves alter- 
nately from each side of a cylinder has been observed in the 
water channel when the stream Mach number was only a 
little above the critical value. Such behavior is to be ex- 
pected if on one side of the cylinder the velocities first 
exceed and then fall below the speed of sound as the flow 
closes in and then breaks away from the surface. With 
increase in Mach number, a strong disturbance originates at 
the edge of the wake approximately 1 diameter behind the 
cylinder and extends into the field of flow (figs. 15 (d), 
18 (d), and 26 (d)). This disturbance oscillates with the 
wake, which is still unstable, and in its incipient stages 
alternates from one side of the model to the other. At still 
higher speed, the flow closes in behind the cylinder so that 
the cylinder has the appearance of a streamline body with a 
strong gull-wing-shaped disturbance at its trailing edge 
(figs. 15 (e), 15 (f), 16 (f), 17 (f), 18 (f), 19 (f), 26 (e), and 
26 (f)). These features appear to be essentially the same 
in the water channel as in the air flow. 




Maximum Mach number , 
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Figure 30 Lines of constant Mach number in the flow field about a 6-inch-diameter circular 
cylinder in the 24-inch water channel at maximum stream Mach number. M~ 0.51. 


Figure 29.— Water-channel maximum test Mach number. 
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EQUIPMENT AND FUTURE DEVELOPMENT 

Measurement of depth at the surface of models and at 
fluid boundaries are best made with pressure orifices and 
burettes because of the capillary rise, whereas the micrometer 
method is suited for measurement of field depths; both 
methods agree very well provided the vertical accelerations 
are small. The burette method is best for obtaining average 
pressures in unstable flow. A method of photography has 
been developed that is simple and completely satisfactory. 
The shadowgraphs of the water flow are strikingly similar 
to the schlieren photographs of the air flow. 

The value of the water channel lies mainly in the low cost 
and convenience of operation. Field surveys are simply 
made, and various features of the flow such as turbulence, 
vortices, separation, and shock formation are easily observed 
and photographed. Streamlines about a model are easily 
obtained by inserting streams of dye in the water ahead of 
the model. 

A larger channel especially designed to minimize boundary- 
layer effects and to secure uniform flow in the test section is 
desirable. The channel should be large enough to permit tests 
at Reynolds numbers above the critical value. A larger 
channel would also be advantageous in reducing the adverse 
effects of the boundary layer or vertical accelerations, or both. 

Additional investigations, both theoretical and experi- 
mental, are needed in order to determine the corrections 
necessary to convert quantities obtained from the water flow 
to the values characteristic of the air flow about correspond- 
ing bodies. 

CONCLUSIONS 

An experimental apparatus and technique have been de- 
veloped for the investigation of the analogy between water 
flow with a free surface and two-dimensional compressible 
gas flow (hydraulic analogy); a preliminary investigation has 
been made and the results of an application of the analogy 
have been presented. The following conclusions are indi- 
cated from this work: 


1. The hydraulic analogy provides a very inexpensive and 
convenient means of investigating high-speed two-dimen- 
sional air flow. The flow may be observed and photo^ 
graphed, and surface and field measurements may be easily 
obtained. Reasonably satisfactory agreement was found 
between the water flow and air flow about corresponding 
bodies, although considerable work in both theory and experi- 
ment is needed in order to convert with quantitative accu- 
racy from the water flow to the flow in air. 

2. With a larger channel, difficulties due to vertical accel- 
erations and subcritical Reynolds numbers might be overcome. 

Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., August 21, 1946. 
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Positive directions of axes and angles (forces and moments) are shown by arrows 


Axis 

Force 
(parallel 
to axis) 
symbol 

Moment about axis 

Angle 

Velocities 

Designation 

Sym- 

bol 

Designation 

Sym- 

bol 

Positive 

direction 

Designa- 

tion 

Sym- 

bol 

Linear 
(compo- 
nent along 
axis) 

Angular 

Longitudinal 

X 

X 

Rolling 

s L 

Y >Z 

Roll 

<t> 

u 

V 

Lateral 

Y 

Y 

Pitching 

M 

z — >x 

Pitch 

Q 

v 

r 

n 

Normal 

Z 

z 

Yawing 

N 

X >Y 

Yaw 


w 

f 











Absolute coefficients of moment 

p ^ p ^ p N 

° l ~([bS Um ~qcS hn ~qbS 

(rolling) (pitching) (yawing) 


Angle of set of control surface (relative to neutral 
position), 8. (Indicate surface by proper subscript.) 


4. PROPELLER SYMBOLS 


D Diameter 

p Geometric pitch 

pjD Pitch ratio 

V' Inflow velocity 

V t Slipstream velocity 

T 

T Thrust, absolute coefficient C T — — 

pn 2 D A 

Q Torque, absolute coefficient G q — 


P 

Os 

v 

n 

4 > 


Power, absolute coefficient C P 


P 

pn*D 5 


Speed-power coefficient = 

Efficiency 
Revolutions per second, rps 



Effective helix angle = tan” 



5. NUMERICAL RELATIONS 


1 hp = 76.04 kg-m/s=550 ft-lb/sec 
1 metric horsepower=0.9863 hp 
1 mph=0.4470 mps 
1 mps=2.2369 mph 


1 lb=0.4536 kg 
1 kg =2. 2046 lb 
1 mi= 1,609.35 m=5,280 ft 
1 m = 3.2808 ft 


